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that  the  peak  wavelengths  of  this  spec- 
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1.  Introduction 

In  1989  Kandpal,  Vaishya,  and  Joshi  of  the 
National  Physical  Laboratory  (India)  published  re- 
sults of  experiments  in  which  they  observed  spec- 
tral shifts  caused  by  a simple  optical  system  [1].  In 
their  experiments  light  from  the  exit  aperture  of  an 
integrating  sphere  was  imaged  by  a two  lens  system 
which  had  an  interference  filter  between  the 
lenses.  They  found  that  when  a small  aperture  was 
placed  in  the  image  plane,  the  peak  wavelength  of 
the  spectrum  of  the  light  measured  on-axis  in  the 
far  zone  of  the  aperture  was  shifted  away  from  the 
peak  wavelength  which  occurred  when  no  aperture 
was  used. 

They  interpreted  these  shifts  as  being  “Wolf 
shifts,”  i.e.,  coherence-induced  shifts  of  the  type 
predicted  by  Wolf  [2-3].  Their  explanation  [4]  for 
the  occurrence  of  the  shifts  was  based  on  their  as- 
sertion (which  was  neither  proven  directly  experi- 
mentally nor  justified  theoretically)  that  the 
filter-lens  combination  eliminated  the  wavelength 


dependence  of  the  complex  degree  of  spectral  co- 
herence of  the  light  in  the  image  plane.  If  this  as- 
sertion is  true,  then  a coherence-induced  change  of 
the  spectrum  is  to  be  expected  when  the  light  prop- 
agates from  the  image  plane  to  the  far  zone,  as  was 
shown  experimentally  by  Morris  and  Faklis  [5]  and 
Faklis  and  Morris  [6].  Kandpal  et  al.  [4]  then 
argued  that  the  introduction  of  a circular  aperture 
in  the  image  plane  helped  in  modifying  the  spec- 
trum in  the  far  zone. 

In  this  paper,  we  will  analyze  theoretically  the 
optical  system  used  in  Ref.  [I].  Many  details  of  the 
system  which  did  not  appear  in  Ref.  [I]  were  pro- 
vided to  us  by  the  group  at  NPL  [7].  The  basic 
outline  of  our  paper  is  as  follows.  In  Sec.  2 the 
optical  system  is  described,  and  the  basic  assump- 
tions to  be  made  throughout  the  paper  are  stated. 
In  Sec.  3 an  approximate  form  for  the  cross-spec- 
tral density  of  the  light  in  the  image  plane  is  ob- 
tained, and  it  is  shown  that  the  corresponding 
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complex  degree  of  spectral  coherence  is  not  wave- 
length independent.  In  Sec.  4 the  spectrum  of  the 
light  on-axis  in  the  observation  plane  is  investi- 
gated, and  an  approximate  form  for  it  is  obtained, 
both  for  the  case  in  which  no  aperture  is  used,  and 
for  the  case  in  which  the  small  aperture  is  used.  It 
is  shown  that,  for  each  interference  filter  used,  the 
peak  wavelength  in  the  latter  case  is  shifted  with  re- 
spect to  the  peak  wavelength  in  the  former  case. 
However,  these  shifts  do  not  agree  with  those  ob- 
served by  the  group  at  NPL.  The  shifts  predicted  by 
our  analysis  are  much  smaller,  so  small  as  to  be  un- 
observable to  within  the  accuracy  of  their  experi- 
ments. In  Sec.  5 our  conclusions  are  presented. 

2.  Optical  System,  Notation  and  Assump- 
tions 

The  optical  system  used  by  the  group  at  NPL  is 
pictured  in  Fig.  1.  Just  to  the  left  of  plane  I there  is 
an  integrating  sphere  of  radius  25  cm  which  has 
a 450  W tungsten  halogen  lamp  at  its  center.  Plane 
I contains  the  exit  aperture,  a,  of  the  integrating 
sphere,  o-  is  a circular  aperture  of  radius 
flo  = 0.12  cm,  and  it  is  in  the  front  focal  plane  of  a 
lens,  Li,  whose  focal  length  is/i  = 5 cm  and  radius  is 
fli  = 0.45  cm. 


Plane  I Plane  II  Plane  III 
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Fig.  1.  The  optical  system. 

There  is  an  interference  filter,  IF,  a distance 
Di  = 20  cm  behind  Li,  and  a second  lens,  L2,  a dis- 
tance Dz  = 30  cm  behind  the  interference  filter.  L2  is 
an  achromatic  doublet,  and  its  focal  length  is/2  = 20 
cm.  The  radius  of  L2  and  the  transverse  dimensions 
of  the  interference  filter  are  all  larger  than  Oi.  It  is 
therefore  a straightforward  matter  to  show  that  Li 
is  the  aperture  stop  for  the  system.  In  their  experi- 
ments the  group  at  NPL  used  six  different  interfer- 
ence filters;  the  shortest  peak  transmission 
wavelength  used  was  422  nm  and  the  longest  was 
652  nm. 


The  exit  aperture  of  the  integrating  sphere  is  im- 
aged onto  plane  11.  It  is  a straightforward  matter  to 
show  that  the  magnification,  M,  for  this  imaging  is: 

M=-/2//,=  -4.  (1) 

Two  different  kinds  of  things  were  done  as  regards 
plane  11.  In  one  case  an  aperture  A of  radius  was 
placed  in  plane  II,  and  experiments  were  done  for 
the  values  = 0.012  cm  and  0.50  cm.*  In  the  second 
case  no  aperture  was  used.  For  notational  conve- 
nience, we  will  treat  the  latter  case  as  if  an  aperture 
A of  radius  0^  = 00  were  in  plane  11.  Also,  since  no 
shift  was  observed  in  the  0.50  cm  aperture  case,  we 
will  not  consider  that  case.  Plane  III,  where  the 
measurements  were  made,  is  a distance  z = 100  cm 
from  plane  11. 

The  following  notation  will  be  used  throughout 
this  paper.  The  locations  with  respect  to  the  optical 
axis  of  points  in  the  planes  I,  II,  and  III  will  be 
specified  by,  respectively,  the  two-dimensional  posi- 
tion vectors  p",  p',  and  p.  The  locations  of  points  in 
the  planes  occupied  by  Li  and  L2  will  be  specified 
by,  respectively,  the  two-dimensional  position  vec- 
tors ai  and  az.  The  lengths  of  vectors  will  be  de- 
noted by  the  corresponding  non-boldfaced  symbols, 

e.g.,  p-  = \p% 

The  following  assumptions  will  be  made.  As  con- 
cerns the  properties  of  the  tight  coming  from  the  in- 
tegrating sphere,  it  will  be  assumed  that  the  exit 
aperture  of  the  integrating  sphere  radiates  as  a uni- 
form lambertian  source.  Let  the  spectral  radiance 
of  this  source,  at  angular  frequency  (o,  be  denoted 
by  The  cross-spectral  density,  IT*'*  (p",  p!z,  (o) 

of  the  light  in  plane  I can  be  written  as  [9-11]:^ 

IT*'*  o))^2ttBo((o) 

xjo{k\f^-p'l\)P,(p'l)Po(f^),  (2) 

where  k is  the  wavenumber  of  the  light, 

k — (x)lc  = 2ttIX.  , (3) 

y'o  is  the  spherical  Bessel  function  of  the  first  kind  of 
order  zero. 


‘ In  Ref.  [1]  the  diameter  of  the  smaller  aperture  was  incorrectly 
reported  as  0.24  cm.  The  correct  value  is  0.024  cm  [8]. 

- In  order  for  Eq.  (2)  to  be  appropriate,  the  radius  of  the  source 
must  be  much  larger  than  the  wavelength  of  the  light.  In  the  NPL 
experiments  the  wavelengths  of  interest  are  visible  and  the 
radius  of  the  exit  aperture  of  the  integrating  sphere  is  0.12  cm, 
so  this  condition  is  fulfilled. 
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jo{u)  = 


sin(w ) 


(4) 


and  Po  is  the  pupil  function  for  the  exit  aperture  of 
the  integrating  sphere, 


Fo(p")  = circ(p7ao).  (5) 

In  Eq.  (3)  c is  the  speed  of  light  in  vacuo  and  A is 
the  wavelength  of  the  light.  In  Eq.  (5)  circ  is  the  cir- 
cle function, 

circ(ii)  = 1,  if 

= 0,  if  |m  I > 1.  (6) 

Secondly,  since  Li  is  the  aperture  stop  for  the  sys- 
tem, it  will  be  assumed  that  the  finiteness  of  the 
transverse  sizes  of  Li  and  the  interference  filter  can 
be  neglected.  Thirdly,  it  will  be  assumed  that  the 
transmittance  of  the  interference  filter  may  be  de- 
scribed by  the  Lissberger-Wilcock  model  [12-13], 
and  that  its  behavior  is  such  that  the  approximation 
described  in  Sec.  4.1  is  valid.  In  addition,  both 
lenses  will  be  treated  as  thin  lenses,  and  the  possi- 
bility of  chromatic  aberration  introduced  by  either 
Li  or  L2  will  not  be  taken  into  account.  The  paraxial 
approximation  will  be  used  throughout  our  work. 


3.  Cross-Spectral  Density  of  the  Light  in 
Plane  11 

In  this  section,  the  cross-spectral  density  in  plane 
II  will  be  investigated.  First  the  coherent  impulse 
response  function  for  the  propagation  from  plane  I 
to  plane  II  will  be  determined.  This  result  will  then 
be  used  to  obtain  an  “exact”  (within  the  paraxial 
approximation)  expression  for  the  cross-spectral 
density  incident  upon  plane  II.  Approximations  ap- 
propriate to  the  NPL  experiments  will  then  be  used 
to  obtain  an  approximate  expression  for  the  cross- 
spectral  density  in  plane  11.  The  corresponding 
spectrum  and  complex  degree  of  spectral  coherence 
in  plane  II  will  then  be  calculated  and  discussed. 

3.1  Coherent  Impulse  Response  Function  for  the 
Propagation  from  Plane  I to  Plane  II 

Let  P"  be  a point  in  plane  I located  by  position 
vector  p",  and  let  P'  be  a point  in  plane  II  located 
by  position  vector  p'.  The  coherent  impulse  re- 
sponse function,  h{p',  p",  m),  for  the  propagation 
from  plane  I to  plane  II  is  the  field  at  P'  due  to  a 
monochromatic,  unit  amplitude  point  source  of  an- 
gular frequency  w at  P".  In  order  to  calculate  this 


function  we  must  first  discuss  the  effect  of  the  inter- 
ference filter. 

The  effect  of  the  interference  filter  on  a plane 
wave  is  depicted  in  Fig.  2.  Here  a unit  amplitude, 
monochromatic  plane  wave  of  angular  frequency  w 
propagating  in  the  direction  specified  by  the  unit 
vector  s is  incident  upon  the  filter  at  an  angle  of  in- 
cidence 6.  If,  as  we  have  assumed,  the  finiteness  of 
the  transverse  size  of  the  filter  can  be  neglected,  the 
effect  of  the  interference  filter  is  to  change  the  am- 
plitude of  the  plane  wave  by  the  factor  t(9,  co),  the 
amplitude  transmission  function  of  the  interference 
filter. 


Fig.  2.  The  effect  of  the  interference  filter  on  a plane  wave. 


When  this  effect  is  taken  into  account,  it  can  be 
shown  (see  Appendix  A)  that  in  the  paraxial 
approximation  h (p',  p",  w)  is  the  product  of  r(p7/2, 
(o),  the  amplitude  transmission  function  of  the  in- 
terference filter  evaluated  at  d-p'/fi,  and  hN(p', 
p",  (o),  the  coherent  impulse  response  function 
(paraxial  form)  for  the  case  when  there  is  no  inter- 
ference filter  present,  i.e., 

h(p',  p",  w)  = t(p'lf2,  uy)hs{p\ p",  <o),  (7) 

where 

^n(p',  p",  co)  = C(a))exp{i/c[(/>(p")  + x(p')]} 


X Besinc 


(kui 

„ 1 , 

P — — p 

M 

and 


\2f2 


M ATf 


(8) 

(9) 


‘/>(p")=/i  + (p"V2/0,  (10) 
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;h;(p')=/2  + D + (l-— )(p'V2/2)  (11)  andf  (pi , p’l,  w)  and  E(p'{,  pi,  w)  are  phase  factors: 

Besinc(u)  = 2/.(u)/u.  (12)  £ (oi.  pJ,  <-)  = exp[a(l -2) , (16) 


In  Eq.  (11)  D is  the  distance  from  Li  to  L2,  i.e., 
D =Di  +D2.  In  Eq.  (12)  J\  is  the  Bessel  function  of 
the  first  kind  of  order  one. 

3.2  Cross-Spectral  Density  of  the  Light  Incident 
Upon  Plane  II 

The  cross-spectral  density  of  the  light  incident 
upon  plane  II,  pi,  (o),  is  related  to  the 

cross-spectral  density  in  the  plane  I by  the  expres- 
sion [14] 

IT'"-)(p;,  p2,(o)  = j I h*(p'u  p'l,  co)  h {pi,  pi,  OJ) 


xiv^^>(p'{,  fH,  cu)dVTdVi,  (13) 

where  h (p',  p",  w)  is  the  coherent  impulse  response 
function  for  the  propagation  from  plane  I to  plane 
II.  Upon  substituting  Eqs.  (2)  and  (7)  into  Eq.  (13), 
and  using  Eq.  (8)  in  the  result,  we  find  that 


E(p'l,  f^,  6;)  = exp[/A:(p^'^ 


(17) 


Equation  (14),  with  F(pl,  p2,  co)  given  by  Eq.  (15), 
is  an  exact  (within  the  paraxial  approximation) 
expression  for  the  cross-spectral  density  incident 
upon  plane  II.  However,  the  integral  on  the  right 
hand  side  of  Eq.  (15)  cannot  be  evaluated  analyti- 
cally. Nevertheless,  for  the  NPL  system  certain 
approximations  are  appropriate  which  simplify 
T(pi',  pi,  co)  considerably. 

First,  it  is  shown  in  Appendix  B that  for  the  NPL 
parameter  values  we  may  set  E{p'l,  pn,  (o)~l  in  the 
integrand  on  the  right  hand  side  of  Eq.  (15).  Hence, 
to  a good  approximation 


F (p'l,  pi,  (o)  = C\(o)E{p[,  pi,  CO 


)// 


X Besincf^ 
Ji 


j Besinc 


P2 


M 


P2 


IF'"-)(p,',  pi,  co)  = 27rBo(co)t*{p{lf2,  (o)t {pi ! fj,  co) 


xF{p[,  pi,  (o)  , (14) 


where 


x;o(A:|A^'-p,"|)dVi'd'p2  . (18) 


Next,  it  is  shown  in  Appendix  C that  for  the  NPL 
parameter  values. 


F(pl,  pi,  (o)  = C\(o)E{p[,  pi,  co) 


X I j F (pT,  p'i,  (o)  Besinc 


p'-jiA) 


X Besinc 


(kai 

V/. 


P2-Pl 


) 


I /Besinc  {f  |pr-^p[|) 


X Besinc 

' /I 


„ 1 , 


)jo{k\f^-p'l\)d^p'ld^f^ 


f2 


xdVi'dVi, 


(15)  xcirc(pi7|M|flo)  circ(p2/|M |flo)  . (19) 
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It  then  follows  from  Eqs.  (14),  (18),  and  (19)  that 
the  cross-spectral  density  of  the  light  incident  upon 
plane  II  is  given,  to  a good  approximation,  by  the 
expression 


IE<")(p;,  P2\  C0)=f^  Boi<o)t*(p[/f2,w) 
h 

xt(p2//2,  (o)E{p[,  P2,  ru)Besinc^^  P2-p[|] 

xcirc(pi7|M|flo)circ(p^/|A/|flo).  (20) 


The  spectrum,  y*""*  (p',  (o),  of  the  light  incident 
upon  plane  II  at  the  point  P'  located  by  the  vector 
p'  and  at  the  angular  frequency  w is  obtained  from 
the  corresponding  cross-spectral  density  by  setting 
P\=p2=p'-^  It  therefore  follows  from  Eq.  (20)  that 


5<"-)(p',  co)  = 

^Boioj)\t{p’/f2,  tu)|'circ(p7|Mlflo).  (21) 


1E™(PI,  fn,  co)=f}  Boico)t*(p[lf2,  oj) 

J2 

Xt{p^/f2,  OJ)  E{p[,  pi,  OJ)  Besincf^ 

' J2 

xcirc(pi7fl)  circ(p2/fl),  (23) 

where  a is  the  smaller  of  \M\ao,  and  a^,  i.e., 
a =min(|M|ao,  Oi). 

The  physical  significance  of  the  parameter  a is 
that  it  is  the  radius  of  the  “secondary  source.”  If 
03  < \M\ao  part  of  the  incident  light  is  blocked  by  the 
aperture  A,  and  therefore  the  radius  of  the  sec- 
ondary source  is  03.  However,  if  03  > |A/|ao  the  aper- 
ture does  not  (within  the  approximation  we  are 
using)  block  any  of  the  light,  therefore  |M|ao  is  the 
radius  of  the  secondary  source.  The  distinction  be- 
tween these  two  cases  is  necessary  because  when 
the  small  aperture  is  used  03  < whereas  when 

no  aperture  is  used  fl3>|M|ao. 

The  spectrum,  S^"’(p',  oj),  of  the  light  in  plane  II 
can  be  calculated  by  using  the  result  in  Eq.  (23): 


It  can  be  seen  from  this  expression  that  this  spec- 
trum is  nonuniform  (since  t depends  on  p')  and  oc- 
cupies the  region  p'<\M\ao,  i.e.,  the  region  which 
corresponds  to  the  geometrical  optics  image  of  a, 
the  exit  aperture  of  the  integrating  sphere. 

3.3  Cross-Spectral  Density,  Spectrum,  and  Com- 
plex Degree  of  Spectral  Coherence  in  Plane  II 

The  cross-speetral  density,  lT‘"’(pI,  p^,  a>)  in 
plane  II  (with  the  effect  of  the  aperture  A taken 
into  account)  is  given  by  the  expression 

(p;> 

(pi,  p2,  co)  circ(pi7fl3)circ(p4/fl3).  (22) 


(p',  6;)  = IE<">  (p',p',co), 

= 7^  Bo(a>)|t(p7/2,  ojf  cnc(p'la).  (24) 

We  see  from  this  expression  that  the  secondary 
source  occupies  a circular  domain  of  radius  a,  and 
that  the  spectrum  of  the  light  inside  the  secondary 
source  is  proportional  to  the  product  of  spectral  ra- 
diance in  plane  I,  B()(a)),  and  transmittance  of  the 
interference  filter,  |t(p7/2,  ^<J)p. 

The  complex  degree  of  spectral  coherence, 
P*''TpI>  P2,  of  the  light  in  the  plane  II  is  defined 
by  the  formula  [16]: 


Upon  substituting  the  approximate  Eq.  (20)  into 
Eq.  (22),  we  find  that  the  corresponding  approxi- 
mate expression  for  the  cross-spectral  density  of  the 
light  in  plane  II  is 


p<">(pl,p[,  a;): 


l^*“’(p!.  P2,  oj) 


VS^"\pl,  c)  V5'"'(p2,  (o)' 


(25) 


^ The  spectrum  S is  also  sometimes  referred  to  as  the  spectral 
density  of  the  light.  It  is  the  optical  intensity  per  unit  angular 
frequency  at  the  position  and  angular  frequency  in  question  [15]. 
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Upon  using  Eqs.  (23)  and  (24)  in  Eq.  (25),  we  find 
that  when  p[<a  and  pn<a,  i.e.,  for  pairs  of  points 
located  inside  the  secondary  source. 
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P2>  = 


E(pi,  P2,  (o) 


t*(p\lf2,  co)  t (pilfi,  (o) 

\t(pi/f2,C0)\  \t  {p^/f2,  0))\ 


xBesinc^^  |p2“p!|)-  (26) 

This  result  shows  that  the  group  at  NPL  was  incorrect 
in  stating  that  the  filter-lens  combination  eliminates 
the  wavelength  dependence  of  the  complex  degree  of 
spectral  coherence  of  the  light  in  plane  II.  In  particu- 
lar, since  the  first  three  factors  on  the  right  hand 
side  of  Eq.  (26)  are  each  unimodular,  it  follows  that 
the  modulus  of  (pi,  p^,  w)  is  equal  to  the  abso- 
lute value  of  Besinc(^fli|p^ -pi'|//2)..  Hence,  the  as- 
sertion by  the  group  at  NPL  that  the  wavelength 
independence  of  p*"’(pi',  p^,  w)  causes  the  spectrum 
to  change  upon  propagation  from  plane  II  to  plane 
III  is  not  correct.  There  may  be  a change  of  the 
spectrum  (we  will  investigate  that  in  the  next  sec- 
tion), but  it  is  not  caused  by  a wavelength  indepen- 
dence of  p,*"*(pi',  p^,  0)). 

In  the  next  section  we  will  refer  to  the  effective 
correlation  length,  L(a)),  of  the  light  inside  the  sec- 
ondary source.  The  correlation  length  of  a Besinc 
type  correlation  function  is  typically  taken  as  the 
smallest  separation,  |p2-p[|,  for  which  the  complex 
degree  of  spectral  coherence  takes  on  the  value 
zero.  Since  Besinc  has  its  first  zero  when  its  argu- 
ment is  equal  to  3.832,  it  follows  from  Eqs.  (26)  and 
(3)  that 

L(6;)  = 0.610^  . (27) 


4.  Spectrum  of  the  Light  On-Axis  in 
Plane  III 

Let  P be  a point  in  plane  III  located  by  position 
vector  p.  When  the  Fresnel  approximation  is  used 
to  propagate  the  field  from  plane  II  to  plane  III,  the 
spectrum  of  the  light  at  P is  given  by  the  expression 
[17] 


5<‘'V  <o)  = 

^ jj  (P'u  P^,  co)cxp[ik(fIi^-p'{y2z] 

xexp[-ik(p2-p[)  • p/z]dfp[dfp2,  (28) 


integrations  extend  formally  from  — oo  to  oo  in  each 
variable.'* 

If  the  secondary  source  obeys  the  Leader  condi- 
tion [18] 


kaL  (co) 
z 


«1, 


(29) 


then  we  are  in  the  far  zone  of  the  secondary  source 
and  the  quadratic  exponential  term  in  Eq.  (28)  can 
be  dropped.  We  will  restrict  our  considerations  to 
the  on-axis  observation  point,  p = 0. 

4.1  No  Aperture  in  Plane  II 

In  this  case  a = |M|flo,  and  the  Leader  condition  is 
not  fulfilled  for  the  NPL  parameter  values.  The  on- 
axis  spectrum  is  therefore  given  by  the  expression 

X exp[//c  {pi ^ - p"i  ^)/2z ]d^pi'  d^p2,  (30) 


where  (pi,  p2,  (o)  is  given  by  Eq.  (23),  with 
a =|M|flo. 

In  order  to  investigate  5*'"’(0,  oj),  let  us  first  in- 
troduce average  and  difference  vectors,  p'  and  p',  in 
plane  II: 


'=Hp2  + p[), 

(31a) 

1 

II 

(31b) 

In  terms  of  these  variables,  it  follows  from  Eqs. 
(30),  (23),  (16)  and  (27)  that 

(0,  = G(p',p',  a.) 

xexp{  — ikp'  • p'/d)d^  p'd^  p' , (32) 

where 

G(p',  p',  <x))  = t*{\p'  -5  p'|//2,  (o) 

'^i(\p'  + I p'|//2>  ^a)  Besinc[3.832p7L(a2)] 

xcirc[(p' -5  p')/|M|fl()]circ[(p' +5  p')/|A/|flo]  , (33) 


where  IT'“’(pI,  p2,  a>)  is  given  by  Eq.  (23),  and  the 


■*  In  practice,  the  circ  functions  in  Eq.  (23)  limit  the  domains  of 
integration  to  the  area  occupied  by  the  secondary  source. 
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and 


1 

d 


1 

z ' 


(34) 


Let  us  now  examine  the  behavior  of  G(p',  p , (o) 
as  a function  of  the  difference  vector  p'  for  a fixed 
value  of  the  average  vector  p'.  The  region  over 
which  the  Besinc  function  is  significant  is  the  do- 
main p'<L((o).  As  p'explores  this  region  the  first 
arguments  oft*  and  t in  Eq.  (33)  each  change  by  at 
most  the  angle 


= 0.610  - . (35) 

J2  Oi 

For  visible  light  this  is  an  extremely  small  angle.  For 
example,  for  A =550  nm  and  ai  = 0.45  cm,  L{a>)l 
/2  = 75  |xrad.  Therefore  we  will  assume  that  the  be- 
havior of  the  interference  filter  is  such  that  when 
evaluating  the  integral  on  the  right  hand  side  of  Eq. 
(32)  we  may,  to  a good  approximation,  neglect  the 
p'  dependence  of  both  t*  and  t,  i.e.,  that  we  may 
approximate  G as 

G(p',  p',  (t))  = \t{p'/f2,  (y)p  Besinc [3. 832p7L (oj)] 


xcirc[(p'  -5p')/|M|uo]circ[(p'  +{p')/\M\ai)].  (36) 


Furthermore,  for  the  NPL  system  parameter  values 
and  the  wavelengths  used  in  their  experiments, 
L{(o)«\M\ao.  Therefore  we  will  now  use  the 
quasihomogeneous  approximation  [19]  which  corre- 
sponds to  replacing  Eq.  (36)  by  the  expression 


= I k(p'//2,  w)p  circ(p7|M|flo) 


X circ(p7b)d^  p'  , 


(39) 


and  b = a id  If 2.  For  the  NPL  parameter  values, 
|M|flo  = 0.48  cm  and  6=0.35  cm.  Therefore,  Eq. 
(39)  simplifies  to 


M{(o) 


w)p  circ(p76)d^  p'. 


2^ 


f|r(0,  cu)P 

-'0 


ede. 


(40) 


where  /3  -blf2. 

Let  us  now  discuss  the  transmittance,  \t(6,  (o)\^, 
of  the  interference  filter.  We  will  assume  that  it  can 
be  described  by  the  Lissberger-Wilcock  model  [12- 
13].  This  model  works  well  for  Fabry-Perot  and  all- 
dielectric interference  filters  for  angles  of  incidence 
which  are  less  than  20°.  It  can  also  be  used  to  de- 
scribe, less  accurately,  double  halfwave  and  induced 
transmission  filters  [20]. 

Let  Ao  denote  the  peak  wavelength  transmitted  by 
the  interference  filter  at  normal  incidence  (0=0), 
(AA)o  denote  its  bandwidth  (FWHM)  at  normal  in- 
cidence, and  To  denote  its  maximum  transmittance. 
According  to  the  Lissberger-Wilcock  model  the 
transmittance  of  the  interference  filter  is  given  by 
the  expression 


G(p',  p',  £U)  = 


t 


(6,  (o)\^  = To 


1 + 


2(A  — Aq) 
. (AA)o 


A(,  ^1' 

(AA)o  77\ 


-1 


|r(p7/2,  w)P  Besinc[3.832p7L(w)]  circ(p7|M|flo). 


(41) 


(37) 

Upon  substituting  the  approximate  form  Eq.  (37) 
for  G into  Eq.  (32),  it  can  be  shown,  after  some 
straightforward  calculations,  that  the  on-axis  spec- 
trum in  plane  III  is  given  by  the  expression 


5<"')(0,  w)  = 


Bo{(o)M{(i>)  , 


(38) 


where  rj  is  the  effective  index  of  refraction  of  the 
interference  filter. 

The  integral  in  Eq.  (40)  can  be  evaluated  (see 
Ref.  [12])  without  any  approximation,  for  the  trans- 
mittance function  Eq.  (41).  The  result  is  that 


M{a)) 


h 

2X 


tan  ' 


2X  1 

.l-l-^(^  + 2Ar)J  ’ 


(42) 


where 


where  A = Ao/37[2ti^(AA)o],  and  ^ = 2(A  - A())/(AA)(). 
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Some  comments  about  the  interference  filter 
parameters  Ao,  (AA)o,  and  t)  are  now  in  order.  The 
effective  indices  of  refraction,  i.e.,  the  values  of  t|, 
for  the  filters  used  in  the  NPL  experiments  are  not 
known.  However,  typical  values  of  t)  are  in  the 
range  1.4<t|<3.4,  with  the  lower  values  being 
more  typical  [21-25].  We  will  therefore  use  -ri  = 2 in 
our  calculations.  The  precise  values  of  Ao  and  (AA)o 
are  also  not  known.  For  reasons  which  will  be  ex- 
plained in  Sec.  4.3,  we  will  take  the  values  of  Ao  for 
the  sb(  filters  to  be  421.9  nm,  484.0  nm,  512.4  nm, 
566.0  nm,  609.1  nm,  and  652.0  nm,  with  their  band- 
widths  being,  respectively,  9 nm,  9 nm,  5 nm,  13  nm, 
8 nm,  and  8 nm. 

4.2  Small  Aperture  in  Plane  II 

In  this  case  03  = 0.012  cm  and  the  Leader  condi- 
tion is  fulfilled.  It  therefore  follows  from  Eqs.  (28) 
and  (23)  that  the  spectrum  can  be  written  as 

A<"'>(0,  w)  = 

j J o,)aiiih. 

A A 

x£(p;,  p^,  OJ)  Besinc(^  |p^ -plljd^i'd'p^,  (43) 

where  A is  the  circular  aperture  of  radius  ax  For  the 
same  reasons  as  were  discussed  in  the  Sec.  4.1,  we 
will  assume  that  the  product  oft*  and  t on  the  right 
hand  side  of  Eq.  (43)  can  be  replaced  by  |/p  evalu- 
ated at  the  average  position,  i.e.,  that  to  a good  ap- 
proximation, 

5<'">  (0,  0))  = 

5o(w)J  j |t(|p;  +P^|/2/2,  w)P  E(p[,  fh,  (v) 

xBesinc(^  |p^-p;|)Tpi'd'p^,  (44) 

For  the  NPL  parameter  values,  it  can  be  shown 
from  Eq.  (41)  that,  to  a very  good  approximation, 
k(|P2  + p[  1/2/2,  w)p«|t(0,  w)p,  for  all  p'l  and  p] 
which  are  in  A.  The  on-axis  spectrum  in  plane  III  is 
therefore  given,  to  a good  approximation,  by  the  ex- 
pression 


5<"''(0,  w)=  ^Bo(co)lt(0,  co)l^N(aj),  (45) 

where 

A A 

xBesinc(^  |p^-pi|)  d'pi'd'p]  . (46) 

This  integral  can  be  evaluated  by  substituting  the 
Fourier  integral  representation  of  the  Besinc  func- 
tion [see  Eq.  (88)]  into  the  integral,  interchanging 
the  orders  of  integration,  and  then  recognizing  the 
resultant  integrals  as  familiar  diffraction  integrals. 
The  result  is  that 


iV(-)=i-i^(^)  Q.C’), 

5=0 

(47) 

where 

(48) 

kuiUi 

f2  ’ 

(49) 

and  the  Q2s{v)  functions  are  those  introduced  by 
Wolf  [26]  and  simplified  by  Petersen  [27]; 

Quiv)=Jl{v)+J](v)  , 

(50) 

02.(1')=-  2^, Vl)  (''[.'4v)J2..,(v) 

-Ji(v)y2i+2(v)]  + 2yyi(v)y2,  + i(v)}  , if  5 >0.  (51) 

4.3  Numerical  Investigation  of  the  Spectrum  On- 
Axis  in  Plane  III;  Comparison  to  the  NPL 
Experimental  Values 

Let  A)  be  the  spectrum,  as  a function  of 

wavelength,  at  the  on-axis  position  in  plane  III,  for 
the  case  in  which  there  is  no  aperture  in  plane  II. 
Let  the  wavelength  at  which  this  spectrum  peaks  be 
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denoted  by  Ap  and  the  bandwidth  (FWHM)  of  this 
spectrum  be  denoted  by  AA.  Let  A)  be  the 

spectrum,  as  a function  of  wavelength,  at  the  on- 
axis  position  in  plane  III,  for  the  case  where  the 
small  aperture  is  in  plane  II.  Let  Ap  denote  the  peak 
wavelength  for  this  spectrum.  The  group  at  NPL 
measured  Ap,  AA,  and  A?  for  six  different  interfer- 
ence filters  [1].  Their  results  are  shown  in  Table  1, 
along  with  the  shift. 


5r(0,  A)=^bo(A)mA(A),  (55) 


where 


Wa(A)  = ToN(2ttcI\) 


1 + 


2(A  — Ap)  ' 
. (AA)o  . J 


(56) 


8Ap  — Ap  — Ap  , (52) 

which  occurred  due  to  the  insertion  of  the  aperture. 


Table  1.  NPL  experimental  results 


Ap,  in  nm 

AA,  in  nm 

Ap,  in  nm 

8Ap,  in  nm 

422.0 

9 

421.0 

-1.0 

484.1 

9 

483.6 

-0.5 

512.4 

5 

514.1 

1.7 

566.1 

13 

564.1 

-2.0 

609.1 

8 

610.3 

1.2 

652.0 

8 

653.2 

1.2 

There  are  several  things  which  should  be  noticed 
about  the  results  shown  in  Table  1.  First,  the  shifts 
are,  in  absolute  value,  of  the  order  of  0.5  nm  to 
2.0  nm.  Secondly,  some  of  the  shifts  are  blueshifts 
and  some  are  redshifts.  Furthermore,  there  are  no 
obvious  trends  in  the  behaviors  of  the  shifts  as  func- 
tions of  either  Ap  or  AA. 

Let  us  now  investigate  the  shifts  that  our  theory 
predicts.  It  follows  from  Eqs.  (38)  and  (42)  that  the 
spectrum  when  there  is  no  aperture  in  plane  II  is 
given  by  the  expression 

5r(0,  A)=^bo(A)mN(A),  (53) 

where  6o(A)  is  the  spectral  radiance,  as  a function  of 
wavelength,  at  the  exit  aperture  of  the  integrating 
sphere  and 


In  order  to  simulate  the  experimental  conditions, 
we  did  the  following  things.  The  color  temperature 
of  the  lamp  was  3200  K [7].  We  therefore  took  the 
spectral  radiance,  6o(A),  to  be  a Planck  distribution 
[28]: 

with  a temperature  of  T = 3200  K.  In  Eq.  (57),  h is 
Planck’s  constant  and  ks  is  Boltzmann’s  constant. 
Furthermore,  we  chose  the  center  frequencies,  Ao, 
and  the  bandwidths,  (AA)o , of  the  interference  fil- 
ters to  be  such  that  the  Ap  and  AA  values  obtained 
for  the  spectrum  5|]"*(0,A)  [from  Eq.  (53)]  agreed 
with  the  experimental  values  of  Table  1 to  the  num- 
ber of  decimal  places  being  displayed  there. 

Table  2 lists  the  interference  filter  parameters, 
and  A<j  and  (AA)(),  we  used  in  our  calculations  and 
gives  the  shifts,  8Ap,  predicted  by  our  theory.  These 
shifts  were  obtained  by  using  a search  routine  to 
find  Ap  and  Ap  from,  respectively,  Eqs.  (55)  and  (53), 
and  then  subtracting  the  two  values. 


Table  2.  Theoretical  values  for  the  shift 


Aji,  in  nm 

(AA)(),  in  nm 

8Ap,  in  nm 

421.9 

9 

0.007 

484.0 

9 

0.009 

512.4 

5 

0.009 

566.0 

13 

0.010 

609.1 

8 

0.011 

652.0 

8 

0.012 

mN(A)  = ^tan  ' 


2X 


l + a^  + 2X) 


(54) 


Here  ^and2f  are  given  by  the  expressions  which  fol- 
low Eq.  (42).  It  follows  from  Eqs.  (45),  (47),  and 
(41)  that  the  spectrum  when  the  small  aperture  is  in 
plane  II  is  given  by  the  expression 


These  shifts  disagree  with  those  measured  by  the 
group  at  NPL.  The  theoretical  shifts  are  approxi- 
mately two  orders  of  magnitude  smaller  than  the  ex- 
perimental shifts.  Furthermore,  the  theoretical 
shifts  are  always  redshifts,  whereas  the  experimen- 
tal shifts  are  in  some  cases  redshifts  and  in  other 
cases  blueshifts. 
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At  this  point  some  further  comments  about  the 
shifts  predicted  by  our  theory  are  in  order.  First,  let 
us  consider  Ap.  A)  is  the  product  of  two  wave- 

length dependent  factors,  /tjaIA)  and  &u(A).  Due  to 
the  presence  of  the  factor  N {lire Ik),  njA(k)  is  not 
centered  at  Ao,  it  is  blueshifted  from  it  by  an  amount 
of  the  order  of  0.0001  nm.  Multiplication  by  the 
Planck  spectrum  then  causes  mA(A ) to  be  redshifted 
by  about  0.1  nm.  Now,  let  us  consider  Ap.  5i^"’(0,  A) 
is  the  product  of  two  wavelength  dependent  factors, 
mN(A)  and  bo(k).  is  not  centered  at  Ao,  it  is 

blueshifted  from  this  value  by  an  amount  of  the  or- 
der of  0.01  nm.  Multiplication  by  the  Planck  spec- 
trum then  causes  mN(A)  to  be  redshifted  by  about 
0.1  nm.  When  we  subtract  Ap  from  Ap  the  two  Planck 
shifts  cancel.  Since  the  blueshift  of  mN(A)  is  much 
larger  than  the  blueshift  of  mA(A),  the  resulting 
shift,  8Ap,  is  a redshift  of  the  order  of  0.01  nm. 

The  physical  origin  of  this  redshift  is  as  follows. 
The  field  incident  upon  the  interference  filter  can 
be  represented  as  a superposition  of  polychromatic 
plane  waves  traveling  in  different  directions,  and 
the  wavelength  at  which  the  transmittance  of  the  in- 
terference filter  peaks  decreases  as  the  angle  of  in- 
cidence of  a plane  wave  increases  [see  Eq.  (41)]. 
Therefore,  for  each  polychromatic  plane  wave  inci- 
dent upon  the  interference  filter  at  an  oblique  angle 
of  incidence,  the  spectrum  of  the  transmitted  light 
has  a peak  wavelength  which  is  less  than  it  would  be 
for  a normally  incident  plane  wave.  In  the  case  in 
which  no  aperture  was  used,  the  contributions  from 
a significant  set  of  such  waves  arrive  at  the  on-axis 
observation  point  in  plane  III.  However,  in  the  case 
in  which  the  small  aperture  was  used,  the  contribu- 
tions from  a large  subset  of  these  waves  were 
blocked  at  plane  II  and  did  not  arrive  at  plane  III. 
As  a result,  the  spectrum  in  plane  III  when  the 
small  aperture  is  used  peaks  at  a longer  wavelength 
than  it  does  when  no  aperture  is  used. 

5.  Conclusion 

Our  paper  contains  two  separate  sets  of  results. 
(1)  In  Sec.  3 an  approximate  form  for  the  cross- 
spectral  density  of  the  light  in  plane  II  was  obtained 
[Eq.  (23)],  and  it  was  shown  that  the  corresponding 
complex  degree  of  spectral  coherence  contradicts 
the  explanation  for  the  shifts  given  in  Ref.  [1].  (2)  In 
Sec.  4 further  approximations  were  made  to  propa- 
gate the  cross-spectral  density  obtained  in  Sec.  3 
from  plane  II  to  plane  III,  and  an  approximate  ex- 
pression for  the  on-axis  spectrum  was  obtained, 
both  for  the  case  in  which  no  aperture  is  used,  and 
for  the  case  in  which  the  small  aperture  is  used.  It 


was  found  that  the  peak  wavelength  of  the  spectrum 
in  the  latter  case  is  shifted  with  respect  to  the  peak 
wavelength  for  the  former  case.  However  the  shifts 
predicted  by  our  analysis  are  much  smaller  than 
those  reported  in  Ref.  [1],  so  small  as  to  be  unob- 
servable to  within  the  accuracy  of  their  experi- 
ments. 

This  brings  us  to  an  important  question.  Our 
analysis  is  predicting  no  observable  shift  (to  within 
the  accuracy  of  the  measurements  of  the  group  at 
NPL),  and  yet  the  group  at  NPL  observed  shifts;  so 
where  is  this  shift  coming  from?  In  our  opinion 
there  are  three  possibilities.  One  possibility  is  that 
their  interference  filters  did  not  behave  in  the  man- 
ner we  have  assumed  in  our  calculations,  and  that 
shifts  arose  as  a result.  A second  possibility  has  to 
do  with  the  spatial  coherence  properties  of  the  light 
in  plane  III.  Since  the  secondary  source  created  in 
plane  II  when  no  aperture  is  used  has  a radius 
which  is  forty  times  as  large  as  in  the  case  in  which 
the  small  aperture  is  used,  the  spatial  coherence 
properties  of  the  light  in  plane  III  will  be  quite  dif- 
ferent in  the  two  cases.  If  this  difference  is  signifi- 
cant enough,  it  may  be  that  the  monochromator 
used  in  the  detection  process  responds  differently 
in  the  two  cases.  A third  possibility,  and  in  our  opin- 
ion the  best  one,  is  that  the  shift  is  caused  by  chro- 
matic aberration  introduced  by  Li,  or  L2,  or  both.  Li 
has  a short  focal  length  and  is  not  an  achromatic 
doublet,  and  is  therefore  an  obvious  candidate  for 
introducing  significant  chromatic  aberration.  L2  is 
an  achromatic  doublet.  However,  it  can  only  be  per- 
fectly achromatic  for  two  wavelengths,  whereas  fil- 
ters with  six  quite  different  center  wavelengths  were 
used  in  the  NPL  experiments. 

Finally,  let  us  comment  on  the  relevance  of  these 
possibilities  for  typical  spectral  irradiance  measure- 
ment systems.  Typical  spectral  irradiance  measure- 
ment systems  do  not  use  interference  filters.  There- 
fore, if  the  shift  is  caused  by  a nonideal  behavior  of 
the  interference  filter,  it  is  irrelevant  for  typical 
spectral  irradiance  measurements.  As  concerns  the 
second  possibility,  it  has  been  known  for  some  time 
that  the  spatial  coherence  properties  of  the  light  in- 
cident upon  a monochromator  do  effect  its  re- 
sponse, viz.,  its  slit  scattering  function  [29-30]. 
Hence  any  comparison  of  experimentally  measured 
spectra  which  have  significantly  different  spatial  co- 
herence properties  at  the  entrance  aperture  of  the 
monochromator  should  take  this  into  account.  Fi- 
nally, if  the  shift  is  caused  by  chromatic  aberration, 
it  will  not  occur  in  spectral  irradiance  systems  which 
use  only  mirrors,  and  is  an  effect  which  should  obvi- 
ously be  taken  into  account  in  any  system  which 
uses  lenses. 
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6.  Appendix  A:  Coherent  Impulse  Re- 
sponse Function  for  the  Propagation 
from  Plane  I to  Plane  II 

6.1  Field  Exiting  Li 

Let  O be  a point  in  the  plane  occupied  by  Li,  and 
let  ai  be  the  two-dimensional  position  vector  in  the 
plane  which  locates  Q.  The  field,  (o),  inci- 

dent upon  Li  at  the  point  Q due  to  a monochro- 
matic, unit  amplitude,  point  source  at  position  P"  in 
plane  I is  given  by  the  expression 


evaluated  at  K =ksi.  It  follows  from  Eqs.  (61)  and 
(63)  that 

co)=  A[^(Sx-^0l)]  , 

lAJl 

= exp[ik(f){p")] 

X Besinc[A:fli|5i -5ui|]  , (64) 

where  soi  == 


U^'-\au  <o)  = 


J_  exp(/A:|ri  -p"|) 
/A  |ri-p"| 


(58) 


6.2  Field  Incident  Upon  La 


where  n = (ai,  /i).  Since  the  angles  involved  are 
small,  this  expression  may  be  approximated 
(paraxial  approximation)  as 

t/<‘-)(a„  a>)=^^^2^exp[/^(a.-p'r/2/,].  (59) 

l Aji 

The  field  exiting  the  lens,  oj),  is  given  by 

the  expression  [31] 


Let  P be  a point  in-between  Li  and  the  interfer- 
ence filter.  Let  r be  a position  vector  which  goes 
from  the  center  of  Li  to  P.  The  field  at  P can  be 
represented  as  a superposition  (angular  spectrum) 
of  plane  waves  [32]: 

L/(r,  oj)  — f a(si,  (o)exp(iks  •r)d^s±.  (65) 

where  s =(sj.,  Sz)  is  a unit  vector. 


[/‘‘■^’(ai,  oj)  = 


Sz  = Vl  -5i  , if 


f/<'-)(ai,  a))exp(-/A:c^/2/i)Pi(ai)  , (60) 


=/V5i-l , if^i>l, 


where  P(a[)  is  the  pupil  function  for  Li, 
Fi(ai)  = circ(ai/fli).  Using  Eq.  (59)  in  Eq.  (60)  yields 


<o)  = 


and 


a{si,  ‘^)=p  U^^^\ksi,  co)  . 


(66) 


(67) 


gxp[/^ip:'l]  exp( •p7/i)/’i(a,)  , (61) 

lAfi 

where 

<A(p")=/.  + (p"V2/i).  (62) 

In  the  next  subsection  we  will  need  &>), 

the  two-dimensional  spatial  Fourier  transform  of 

<o), 

w)  = 

/ (o)  exp(  — iK  • ai)d^ai,  (63) 


In  Eq.  (66),  Sa.  = |sx|. 

In  spherical  polar  coordinates,  as  measured  from 
the  center  of  Li,  s=(l,  6,  4>).  In  this  notation 
si  = (5i,  5y)  = (sin0cos</),  sinPsint^).  Therefore  the 
angle  of  incidence  of  the  plane  wave  exp(/^s  t) 
with  respect  to  the  normal  to  the  surface  of  the  IF 
is  6 [see  Fig.  2],  where  s±  = sin  6~  6. 

Let  t{6,  co)  be  the  amplitude  transmission  func- 
tion of  the  filter  for  a plane  wave  of  angular  fre- 
quency CO  and  angle  of  incidence  6.  If  we  neglect  the 
finiteness  of  the  transverse  size  of  the  IF,  the  effect 
of  the  interference  filter  is  to  change  the  amplitude 
of  each  plane  wave,  exp(/A:5  -r),  by  the  factor  t{s±, 
co).  The  amplitude,  call  it  d(si,  co),  of  the  plane 
wave  exiting  the  filter  is  therefore 
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d(si,  oj)=  t(s±,  (o)a{si,  0))  , 

= jjtisu  oj)  (o)  , 


hip',  p",  oi)  = 


exp{/^[<A(p")  + A^(p')]} 


= ^ tisu  ^o)  exp[/7:0(p")] 

X Besinc[/:rJi|si -soi|]  . (68) 

Let  r2  = (a2,  L))  be  the  position  vector  which  goes 
from  the  center  of  Li  to  the  point  in  L2  located  by 
the  two-dimensional  position  vector  aj.  The  field, 
U^^~\a2,  (o)  incident  on  L2  can  be  represented  as  a 
superposition  of  plane  waves: 

U^^~\a2,  (o)=  f d(sx,  co)exp(i/cs  ‘r)d^s±, 

= f d(s±,  0))  Qxpiiksi  • a2) 


X Besinc(/:ai|(p7/2)-(p7/i)|)t(p7/2,  w)  , (74) 
where 

A'(p')=/2+D  + (i-^)(p'V2/2).  (75) 

The  magnification  of  the  imaging  from  plane  I to 
plane  II  is  M = -/7/i.  It  therefore  follows  from  Eq. 
(74)  that  the  coherent  impulse  response  function 
can  be  written  as 

hip',  p",  a))  = tip'lf2,  ft>)/jN(p',  p",  (o),  (76) 


X exp(//cszD)d^5'j..  (69) 

Since  the  angles  involved  are  small,  we  may  make 
the  approximation  = Vl  - ~ 1 - 5 in  Eq.  (69), 

and  we  find  that 


where  ht<ip' , p",  (o)  is  the  coherent  impulse  re- 
sponse function  if  there  is  no  interference  filter, 

/in(p',  p",  <o)  = Cico)cxp{ik[<j>ip")  + xip')]} 


f/(2- 

^(0-2,  w)  = exp(i7:D)/u(5i,  co)  expiiks±  '02) 
'X&xpi-ikDs\l2)d‘Si_.  (70) 


X Besinc 


(kai 

V/. 


tt 


P 


ill) 


6.3  Final  Expression 

The  coherent  impulse  response,  hip',  p",  (o),  for 
the  propagation  from  plane  I to  plane  II  is  the  field 
in  plane  II  due  to  the  point  source  in  plane  I,  i.e., 

hip',p",(o)  = U^yp',co)  . (71) 

If  we  neglect  the  finite  transverse  size  of  L2,  this 
equation  can  be  written  as  [33] 


hip',  p",  U^^-\kp'lf2,  co),  ill) 

IAJ2 


where  tA(p')  =/2-l-  (p'72/2).  It  follows  from  Eq.  (70) 
that 


D^^-\ksu  0J)  = 

A^fl(si,  (o)  txpiikD)  Qxpi  — ikDs\l2).  (73) 


Here  (hip")  and  xip')  are  given  by  Eqs.  (62)  and 
(75)  respectively,  and 


Cio>)  = 


J_  TTfll 

M • 


(78) 


7.  Appendix  B:  Derivation  of 
E ip'ij  fH,  (o)~l 


Let  us  first  define  average  and  difference  vari- 
ables in  plane  I as,  respectively, 

P"=iip'i  + P'l),  (79a) 

p"=Pi-p'l.  (79b) 


In  terms  of  these  variables  Eip'l,  fn,  (o)  [see  Eq. 
(17)]  can  be  rewritten  as 


Upon  substituting  Eq.  (73)  in  Eq.  (72)  and  Eq. 
(68)  in  the  result  we  find  that 


Eip'l,  p{,  w)  = exp 
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As  pi'  and  fin  each  explore  the  domain  a,  the  max- 
imum value  of  p"  = |p"|  is  ao-  In  the  integrand  of  Eq. 
(15)  the  range  of  values  of  p"  = fn—p'{  over  which  jo 
is  significant  is  p"<XI2.  Hence  as  pi'  and  pn  range 
over  the  domain  where  the  integrand  is  significant, 


(81) 


If  Ki  » 1 then  the  first  and  second  Besinc  func- 
tions in  the  integrand  on  the  right  hand  side  of  Eq. 
(85)  are,  respectively,  very  sharply  peaked  about  the 
values  ui=vi  = p'JMao  and  U2  = V2  = pilMuo.  As  a re- 
sult, if  vi  or  V2  lies  outside  the  unit  disk,  then  the 
integral  on  the  right  hand  side  of  Eq.  (85)  is  approx- 
imately zero,  and  if  vi  and  V2  both  lie  within  the  unit 
disk  the  integrals  in  Eq.  (85)  may  be  approximated 
by  extending  the  limits  to  ± oo,  i.e.. 


For  the  NPL  system  flo  = 0.12  cm  and  /i  = 5 cm; 
therefore  over  the  domain  where  the  integrand  is 
significant  the  real  and  imaginary  parts  of  E{pi,  fn, 
(o)  fulfill,  respectively,  the  inequalities: 


B(p'i,  p^,  (o)~ao  circ(pi'/A/flo)  c\xc(p2/Mao) 
x/J  Besinc(Ki|Mi -vi|)  Besinc(Ki|M2  — V2I) 


0.997  < Re[£(p'i,  ^',  w)]  < 1 , (82a) 


xyo(K|M2-«i|)d^Wid^W2  . (87) 


-0.075  < Im[£(p'i',  pi',  w)]  < 0.075  (82b) 


8.  Appendix  C:  Derivation  of  Eq.  (19) 

Let 


The  right  hand  side  of  Eq.  (87)  may  be  simplified 
by  using  the  Fourier  representations  of  the  func- 
tions in  the  integrand: 


Besinc(/<:i|«;  - Vj  |)  = 


B(pi,  pi.  a,)  = |jBesinc(^ 


X Besinc 


(ka 
.nc[j 


/. 

)7o(^|p::'-p'.'|)  dV'i'd^pi, 

(83) 

and  let  us  make  the  change  of  variables 

Uj  = p"/ao,  (j  = l,2)  (84a) 

Vj  = pj'IMao,  (J  = 1,2)  (84b) 

In  terms  of  these  new  variables,  Eq.  (83)  can  be 
rewritten  as 

B{p[,  p2,  a))  = ao  j j Besinc(Ki|Mi -vi|) 

u.d.  u.d. 

xBesinc(K’i|u2-V2|)/o(x|«2-Mi|)d^uid^W2  , (85) 
where  u.d.  is  the  unit  disk,  u.d.  -{u:u<l},  and 

(86a) 
(86b) 


Kl  = 


koiUo 


h ’ 

K = kao  . 


X exp[  - iKj  • (uj  - V; )]  d^  /C; , (/■  = 1 , 2),  (88) 

yo(K|M2-M||)  = 


If  drc{KlK)  exp[-/A:  • (u2-Mi)]d^/C 

27rK-J  V1-(/C/k)2 


(89) 


Upon  substituting  Eqs.  (88)  and  (89)  into  Eq.  (87), 
we  find,  after  performing  some  simple  integrations, 
that 

B{p[,  p2,  oj)=D(io)  circ(pi'/A7flo)  c\xc{fhJMao) 

X f circ(A:/x.)  circ(A:/K) 

= D{(o)  circ(pI/Mflo)  circ(p2/A/flo)  f 


where 


X exp  [iK  ' (v2 - vi)]d^K  , 


, (IngoY 


(90) 


(91) 
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and  the  step  performed  in  Eq.  (90)  follows  from  the 
fact  that  ki/k  =fli//i  < 1. 

If  (fli//i)^<<  1,  Eq.  (90)  can  be  approximated  as 
B{p[,  p2,  (x))^D{(o)  circ(pi7Mflo)  circ(p^/Mflu) 

xj  ciTc(K/Ki)exp[iK  ’ {v2  — vi)d^K, 

= D((o)  circ(pi7Mflo)  circ(p2/Afao)77K7 

X Besinc(Ki|v2-Vi|), 

"2CYJ)  Besinc(^Ip^-p;|)  circ(p;/|M|ao) 

xcirc(p^/|A/lflo)  . (92) 

For  the  NPL  system,  flo  = 0.12  cm  and  fli//i  = 0.09, 
and  the  peak  wavelengths  of  the  interference 
filters  used  in  their  experiments  were  in  the  range 
422  nm  to  652  nm.  For  these  parameter 

values,  1.04  X 10^  < < 1.61  x 10^  and  (fli//i)^ 

Ji 

= 8.10x  10“^;  therefore,  ki>>  1 and  {a\lf\y<<  1.  It 
then  follows  that  the  approximation  Eq.  (92)  is  ap- 
propriate for  the  NPL  system. 
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